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1.0 Introduction 


This report summarizes work undertaken by Intermetrics In support of 
the JPl Mars Sample Return (MSR) mission feasibility evaluation. The analysis 
has concerned itself with evaluaMng alternative methods of rec^/erlng a sample 
module from a trans-earth trajectory originating in the vicinity of Ma,‘S. 

The major modes evaluated are.* 

(1) Direct atmospheric entry from trans-earth trajectory. 

(ii) Earth orbit insertion by retro- propul si on. 

(iii) Atmospheric braking to a capture orbit. 

In addition, the question of guided vs. unguided entry vehicles has been 
considered, as well as alternate methods of recovery after orbit insertion 
for modes (ii) and (iii). 

Chapter 2 presents a summary of results and conclusions reached in sub- 
sequent chapters. Chapters 3 and 4 discuss analytical results for aerodynamic 
and propulsive maneuvering vehicles respectively. Chapter 5 discusses system 
performance requirements, and Chapter 6 alternatives for inertial systems 
implementation. Chapter 8 discusses orbital recovery operationsand Chapter 
9 describes further studies required to resolve the recovery mode issue. 
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2,0 Summary and Conclusions 


Based on the data presented in Chapters 5 and 6, all of the candidate 
methods appear to be technically feasible. The most severe system performance 
requirements are presented by low-angle, guided entry vehicles, and they apply 
primarily to the inertial subsystem initial alignment. Failure to achieve 
the required alignment accuracy, or significant degradation of the inertial 
subcomponents will cause large dispersions in impact point, and possibly 
loss of the sample module. Because an unguided vehicle does not depend on 
the correct operation of guidance, navigation, and control subsystems, it is 
not vulnerable to the class of failure modes arising out of the presence of 
such systems. For this reason, unguided aerodynamic vehicles are considered 
more reliable than guided, provided they can achieve the mission objectives. 

One such unguided vehicle is the high flight path angle, unguided ballistic 
entry body, which is the recommended approach to the minimum sample return mission. 

Any vehicle targeted to enter or fly through any portion of the terrestrial 
atmosphere is vulnerable to on-board or ground system failures which prevent 
attainment of the required trajectory conditions or vehicle flight regime. 

If such failures occur prior to the time when the vehicle is first committed 
to an atmosphere encounter trajectory, the only result is loss of the sample 
module into Interplanetary space. After committment to such an encounter 
trajectory, however, such failures must be considered as equivalent to 
compromise of the biological integrity protocol associated with the recovery 
phase. Therefore the optimum recovery mode is considered to be one which does 
not require such tarqeting. This implies an earth orbit insetion by retro- 
propulsion followed by an orbital recovery system. It is suggested here that 
remote earth orbital operations, such as retrieval of objects above Space 
Shuttle operational altitudes, is likely to be a requirement shared by a community 
of users during the 1980's to 1990's time frame, and that development of 
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Shuttle compatible equipment should be undertaken. It is further suggested 
that return of the sample to a terrestrial laboratory can be accomplished by 
use of the Shuttle, provided that special measures are taken to ensure 
integrity of the sample container against an entry catastrophe. 

Further studies are recommended to define structural, thermal protection, 
and weight characteristics of the "minimum" sample return module; to evaluate 
pre-deployment failure and trajectory control problems for the minimum mission; 
to develop earth encounter constraints for safe earth orbit Injections for the 
optimiun mission; and to assess the size and requirements of the community 
interested in remote earth orbital operations. 
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3.0 Aerodynamic Analysis 


3.1 Introduction 

Earth based retrieval of the Mars sample requires dissipation of the high 
kinetic energy acquired by the vehicle on the return orbit. A feasible 
approach to effecting this deceleration lies in exploiting the aerodynamic 
forces induced by motion of a vehicle through the Earth's atmosphere. 

When the vehicle enters the atmosphere, it will experience drag, acting 
in direction opposite to the vehicle's motion, and possibly lift, acting 
normal to the direction of motion. These two forces, in combination with 
the Earth’s gravitational attraction, centrifugal force, and the vehicle's 
inertia will completely condition the ensuing trajectory. 

A vehicle encountering the atmosphere may sliCw through it and re-exit 
further along the path, or it may fall to the surface in a path wholly 
contained within the atmosphere. The occurrence of one or the other will 
depend on the vehicle's altitude rate at atmospheric encounter and the time 
history of altitude acceleration as governed by the radial component of the 
resultant of the above four operating forces on the vehicle. 

For the re-exit trajectories, a critical consideration is the magnitude 
of the re-exit velocity. If the vehicle exits at hyperbolic speed, it will be 
lost into interplanetary space; at sub-parabolic exit speeds, the vehicle 
will be captured in an elliptic orbit around the earth. !;, the latter case 
the vehicle will re-enter the atmosphere in following orbital revolutions, 
suffering further deceleration and eventually falling to a landing in the 

earth's surface. 

The above trajectory characteristics are defined by the vehicle s aero- 
dynamic coefficients and the dynamic conditions at atmospheric encounter. For 
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given vehicle aerodynamic capabilities, ranges of initial conditions (entry 
corridors) can be defined which will yield corresponding terminal conditions: 
hyperbolic re-exit, re-exit to captured elliptical orbit (atmospheric braking), 
or direct fall to earth landing (direct entry). 

When the returning vehicle encounters the atmosphere at orbital speeds, 
a shock wave will form ahead of the vehicle's nose, heating the local 
atmospheric section to very high temperatures. As the vehicle plunges deeper 
into a denser atmosphere, it will increasingly be heated by this enveloping 
layer, while the speed of the vehicle will continuously be reduced by the ^ 
braking action of the atmosphere. In this fashion, a portion of the vehicle s 
kinetic energy is converted into heat. 

The magnitude of the instantaneous deceleration applied to the vehicle, 
and the rate of heat flow to the body represent vehicle design parameters 
which must be met to insure vehicle survival during re-entry. 

The operating drag deceleration varies directly with atmospheric 
density p and the square of vehicle velocity, v 2 . As the path penetrates 
the atmosphere, the density increases rapidly, while velocity decreases due 
to the drag effecc. The combination of these two countervailing effects 
causes a deceleration peak, with subsequent decreasing decelerations. 

Heat rates into the vehicle, assuming a blunt nose design, result 
primarily from convective and radiative transfers from the shock layer. 

The operating rates, analogous to the vehicle deceleration, vary with 
atmospheric density and vehicle velocity, and exhibit a similar peak in 
their time histories. High heat rate trajectories are associated with high 
peak decelerations; however, the total heat flow into the body throughout the 
oath tends to be smaller. That is, the high decelerations cause a rapid fall- 
off in vehicle velocity, and a consequent rapid reduction in the operating heat 
rates. The total heat flow, which reflects mainly this high, but short 
duration peak, is normally lower than for oaths exhibiting lower peaks winch 

subsist for longer time periods. 

An important consideration in the design of any re-entry scheme is the 
sensitivity to errors, both in the nominal atmospheric encounter dynamic 
conditions and in the assumed vehicle and environment parameters, it is 
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particularly so in the case of hyperbolic return orbits because of the 
narrowness of the potential entry corridors, and the disastrous consequences 
of exceeding them. That is, on one side, hyperbolic re-exit and loss of the 
vehicle is possible; on the other, structural failure or -turn-up of the 
vehicle in the atmosphere may result. Even if these extremes are avoided, 
significant dispersions in the final capture conditions may make eventual 
recovery of the sample difficult, if not impossible. 

Tne above topics are analyzed in the following sections. Section 3.2 
discusses the baseline Mars return orbits in terms of their implications for 
the atmospheric encounter problem, and identifies the main control parameters 
utilized in this study. Background material for this and following sections 
’s presented in Appendix A in the form of an analysis of the equations of 
motion for a venicle moving within the earth's atmosphere. 

Se ;ion 3.3 discusses the direct entry capture mode. We first consider 
a ballistic vehicle and evaluate entry corrido*\ ranging capabilities 
expectec vehicle stresses, and error sensitivities. We then consider the 
effect of major modifications to the vehicle's capabilities, in particular, 
the addition of a lifting capability. 

Section 3.4 considers the atmospheric braking capture mode. A similar 
format as the previous section is followed, with discussion of terminel 
conditions, vehicle stresses, and error sensitivities for a ballistic vehicle, 
and extensions to a lifting vehicle. Finally, we briefly discuss the decay 
characteristics of the established capture orbits. 

3 . 2 Atmospheric Encounter Condition s 

For a vehicle possessing a stable configuration at zero roll angle all 
operating forces (i.e., drag, lift, gravity, centrifugal) are contained in the 
plane defined by the vehicle's velocity and position-relative to earth-center 
vectors. The integration of the equations of motion in a wind-relative set 
of coordinate axes (see Appendix A) can then proceed from specification of the 
vehicle's earth relative velocity v and flight path angle 6 at the atmospheric 
encounter altitude h. 
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For the Mars return hyperbolic orbits, the vehicle velocity and associated 
flight path angle, relative to a non-rotating earth, are given In terms of the 
orbit hyperbolic excess velocity v w and the vacuum perigee radius r . The 
encounter azimuth cannot be so freely specified, but rather Is constrained 
by the possible combinations of entry latitude L and orbital Inclination 1. 

The later Is, furthermore, limited from below by the declination of the 
hyperbolic asymptote. However, the range of entry conditions can be investi- 
gated by arbitral’ ly assuming values of 0° » -90°; L » 0°) and 180° 

* 90°; L « 0°) for the Inclination. The equations of motion Involve, however 
atmospheric alrmass relative parameters, rather than Inertial quantities. The 
two variable sets are related through Coriolis effects. 

i ' V (3-1) 

where O Is the earth's angular velocity vector and r is the vehicle position 
vector with respect to the center of the earth. 

Sample Mars return orbits, as described by the hyperbolic excess velocity 
vector v^, are given in Table 3-1. 


Table 3-1: MSR Sample Orbits 



> 


where X Is the right ascension and 6 the declination. The corresponding 
atmospheric encounter velocities are given In Table 3-2. 



i 
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Table 3-2: MSR Atmospheric Encounter Velocities 



v i (kms/sec) 

, max (kms/sec >* 

v m1n (kms/!ec >* 

1 

11.49216 

11.7012 

11.27927 

2 

11.53876 

11.6823 

11.39341 

3 

11.42793 

11.5987 

11.25453 


*Note: for Y. = 0 


The first column represents encounter velorit es relative to a non-rotating earth, 
the last two atmospheric airmass relative velocities for maximum and minimum 
relative speed encounters. 

Thus, expected hyperbolic excess velocity variations do not significantly 
affect the atmospheric encounter speed. On the other hand, velocity dispersions 
at Mars departure are capable of producing significant changes in the orbit 
perigee altitude, and hence, in the atmospheric encounter flight path angle. 

This variable shall therefore constitute the control parameter for the following 
investigations of the atmospheric entry process, and will be utilized as the 
basis for definition of entry corridors for the considered vehicle capture 
modes. 

The altitude at which aerodynamic flight is initiated is defined to be 
that of the outer edge of the earth's feasible atmosphere. This boundary is 
selected to be at an altitude above the earth's surface of 121.92 kms 
(400,000 ft), at which point the atmospheric density is of the order of 
1.87 x 10" 8 kgs/m 3 . 
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3.3 Direct Entry 

We are concerned here with the capture characteristics for direct 
atmospheric entry of unguided, fixed L/D vehicles. As discussed in Section 
3.1, the parameters of concern are listed as 

. entry corridor 
. landing footprint 
. maximum g-load 
. maximum heat rate 
. totai heat load 
. error sensitivities 

Section 3.3.1 evaluates the above items for a ballistic vehicle, parameter- 
ized on the atmospheric encounter flight path angle and velocity. Section 
3.3.2 considers the impact of changes in the vehicle's ballistic coefficient, 
and Section 3.3.3 discusses the effect of vehicle lifting capabilities. 

3.3.1 Ballistic Vehicle 

lhe b.sellne vehicle is e sphere of 1 meter radius, 30 kgs of mass, aod 
ballistic coefficient B of .25 m Z /kg (i.e., Cp = 2.4). 

Figure 3-1 presents plots of down-range in kms versus earth relative 
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DIRECT ENTRY 



FLIGHT PATH ANGLE - DEGREES 


















flight path angle in degrees. Curves are presented both for maximum and 
minimum earth relative velocity entry. 

The plots are seen to rise steeply as the flight path angle approaches 
-4°. This unbounded rise indicates the presence of the overshoot boundary 
for the direct entry corridor - that is, at lower flight path angles, the 
vehicle does not fall to a landing, but rather re-exits from the atmosphere. 

The overshoot boundary location is governed essentially by the vehicle’s 
Initial altitude rate; therefore, capture for maximum relative velocity 
entry is possible at shallower encounter flight path angles than for the 
minimum relative velocity case. Table 3-3 presents sample comparative values 
of achievable down-ranges. 

The landing footprint, i.e., the range of reachable landing latitudes 
and longitudes, is determined by the downrange and the atmospheric encounter 
latitude, longitude and azimuth. 

Due to the rotation of the earth, the encounter longitude is determined 
primarily by the vehicle time of arrival at the encounter altitude. The 
corresponding latitude and azimuth, however, cannot be so freely specified. 

Specification of the return orbit vacuum perigee altitude (or 
equivalently, the encounter inertial flight path angle) fixes the angular 
deviation between the encounter radius and the direction of the hyperbolic 
excess velocity vector v*. Referring to Figure 3-2, the angle v between the 
hyperbola's asymptote and the direction of the perigee radius is given by 

v = cos -1 (l/e) ( 3 . 2 ) 

where the eccentricity e is defined by 


e * 1 + r p /a 


while the polar angle f is given by 
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Figure 3-2: Hyperbolic Encounter Geometry 
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Direct Entry Characteristics 
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(3-4) 


f = cos"^[((r p / r)(l+ e ) - l)/ e ] 

The total angle S - f + v is therefore defined in terns of , and r (or v 1 
and -s Plotted in figure 3-3 versus the encounter inertia, flight pa’t ’ 
angle for the first return orbit listed in Table 3-1 (,.e.. v * * . 

table'do 1 " *7 Vanat '°" S f ° r the rernai " tn9 rt!tu, '" trajectories" Isted In the 

table do not materially change the resulting 9 angle. 

veinrir a " 9le 9 defi " M * C °" e 0f P0sni0n about the hyperbolic excess 
velocity vector direction, with vertex at the earth center. The intersection 

.tmosle C ° ne W ' th ^ Mrth - CCTtered s P he " ° f radius corresponding to the 

rr c e " C0Unter aU,tUde " 4 c1rcle - “ h,ch constitutes the locus of 
atmospheric encounter points for a given flight path angle and hyperbolic 

within r WeCtdr ' add,tl0n ' a ” enc °unter points must be contained 
within a return orbital plane, which also includ-s the v. vector. Based on 

two geometric considerations, and utilizing standard spherical 

tngonometry relations, the entry latitude and inertia, azimuth can be shown 


s1nL e * sini cose + sineysin 2 , - sin 2 S 


sin 0 ^. = -cos i / cos L 


(3-5) 


:;::r ion ° f the vector and 1 «• - «- 

NOte ,r ™ tba ‘hat. for inclination angles sy™tric about 90” 

e obt i'n ' ’ ' atUu<,e bU ‘ a "‘ 1s «tric entry azimuths 

! furthermore, for any orbital inclination, two possible entry 

the sro^rr: correspondins to ° pp ° sue directi ° ns ° f -«« «“h<" 

constrained 2 ha " d ' "'° t,0n ° f the veh1c1e ‘he atmosphere Is 

constrained within an earth-fixed plane. The landing latitude is then given 
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sinL^ * sinL g cos(R/r e ) + cos L g sin(R/r e )cos ^ (3-5) 

where R is the downrange and \p is the earth-relative velocity-based azimuth 
angle. 

The band of reachable landing latitudes is then determined by specifica- 
tion of v B , i, and figure 3-4 presents plots of landing point latitudes 
as a function of flight path angle for various orbital plane inclinations. 

The return trajectory represents the first case listed in Table 3-1, i.e., 
v w 8 3.07 km/sec, <5 « 27.13°. 

In the generation of these curves, the variation of achievable down- 
range with entry azimuth for a given flight path angle, was neglected. As 
shown in Figure 3-1, this effect is small for all encounter flight path 
angles not in the immediate vicinity of the overshoot boundary. In addition, 
the difference between inertial and relative encounter azimuth was also 
neglected. Again, this effect is not significant except in the Immediate 
vicinity of the overshoot boundary. 

Figure 3-5 presents the bounds of reachable latitudes for the two 
remaining return trajectories in Table 3-1, i.e., v w * 3.24 km/sec, 

6 8 -52.45° and v^ 8 2.32 km/sec, 6 8 -43.44°. These bounds were obtained 
by solving for the landing point latitude for the 90° inclined orbit. 

A characteristic of all direct entry paths considered is that they 
consist of a relatively brief period of deceleration followed by a similar 
equilibrium descent. That is, following the peak deceleration point, the 
vehicle is captured in a vertical air column, and descends in a regime where 
the drag deceleration is equal and opposite to gravity. The vehicle's 
velocity at any point is only a function of altitude, and is independent of 
initial encounter conditions. 
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The vehicle terminal velocity can be explicitly solved for in terms of 
s-rface atmospheric density and gravity magnitude, and is 


v f = 7.98 m/sec 


(3-7) 


for the baseline vehicle. Furthermore, since most of the path is comprised 
by this equilibrium fall, the time-of-f light from atmospheric encounter to 
landing is relatively independent of encounter conditions. Sample vc.ues of 
this parameter are also listed in Table 3-3. 

Figure 3-6 presents peak g-loading vs. encounter flight path angle 
curves for maximum and minimum relative velocity entry. The minimum P^k 
loading for direct ballistic entry is seen to be on the order of 5 g’s at 
the overshoot boundary, rising to 325 g‘s for vertical entry. 


The maximum peak loading for maximum relative velocity entry does not 
occur at the earth relative flight path angle of -90°. Rather, it occurs 
at an inertial flight path angle of -90°, which, due to the motion of the 
earth, corresponds to a relative angle of about -87.5°. The relative angle 
of -90° is achieved at a slightly lover inertial angle on a positive heading 
which results in partial cancellation of inertial velocity by earth speed, 
and therefore, a lower velocity of the vehicle relative to the air mass. 


Figure 3-7 presents peak stagnation point heating rates and total heat 
load on the vehicle as a function of encounter flight path angle for maximum 
and minimum relative velocity entry. As shown, peak heating rates increase 
as the flight path angle tends towards -90°; however, as discussed in 
Section 3.1, total head loads on the vehicle decrease. The heating rate 
curve for maximum relative velocity entry exhibits, similarly to the g-loading 
curve, a maximum at an angle offset from -90°, and for the same reasons. 

Table 3-3 includes sample values of the above vehicle stress parameters for 
various encounter flight path angles. 
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Figure 3-8 and Table 3-4 present down-range dispersion sensitivities 
to various error sources. Four driving errors are considered: 

. error in assumed encounter velocity 
. error iri assumed encounter flight path angle 
. error in assumed vehicle ballistic coefficient 
. error in assumed atmospheric density profile 

The provided curves represent the maximum relative velocity entry case. The 
similar performance exhibited in Figures 3-1, 3-6, and 3-7 by the maximum 
and minimum relative velocity entries indicates that similar qualitative 
behavior as shown in Figure 3-8 would be expected for the unrepresented 
cases. 

The errors considered represent biases on the entry parameters. The 
ballistic coefficient error reflects uncertainties as to the vehicle's 
drag coefficient, reference area or mass. The atmospheric density error 
represents uncertainty in the density profile of a vertical air column; the 
formulated error is in the form of a percentage of the seasonal deformation 
AP of this profile, which tends to redistribute the air mass along the column. 

The general shape of all four curves exhibits high sensitivities close 
to the overshoot boundary, which rapidly decay as the encounter flight path 
angle tends to -90°. 

Specification of the most critical error source depends also on the 
definition of the level of expected driving errors. For normally expected 
values, it is clear that the encounter flight path angle represents the 
highest sensitivity error sources. 

An alternative evaluation of these curves would lead to a specification 
of maximum allowable error source uncertainties for a specified level of 
acceptable down-range dispersions. Assuming uncorrelated error sources, the 
down-range error variance is given by 
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Figure 3' 


. Down Range Error Sensitivities - Direct Entry 
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Baseline Ballistic Vehicle Down-Range Error Sensitivities 














o R 2 s (3R/3v) 2 o y 2 + (3R/3y) 2 o^ 2 + (3R/3B) 2 0g 2 + (3R/3p) 2 a p 2 


(3-8) 


Assuming equal weights for all error sources, the acceptable uncertainty 
levels are given by 


o y * a R /2(3R/3v) 
c y * o R /2(3R/3y) 
o B * o r /2(3R/3B) 


(3-9) 


o p a o r /2(3R/3p) 

Given the down-range dispersion level (o R ), these functions can readily 
be generated from the curves of Figure 3-8. 

The figures and tables presented above provide a basis for evaluation 
of the viability of the direct ballistic atmospheric entry mode for capture 
of the returning Mars sample. Selection of acceptable entry trajectories is 
governed by 


. maximum allowable vehicle stresses (g-loading, heating) 

. desired landing footprint 
. acceptable landing dispersions 

Figures 3-6 and 3-7 allow translation of specifications of acceptable 
vehicle stresses into a range of permissible encounter flight path angles, 
and constrains the candidate values towards the horizontal or overshoot 
boundary. On the other hand, smaller flight path angles imply higher error 
sensitivities, and hence, higher landing point dispersions. For a given set 
of encounter condition uncertainties, constraints on acceptable landing 
dispersions will limit the range of permissible encounter flight path angles 
from below (i.e., constrains them towards the vertical), as shown in 
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Figure 3-8. Finally, the desired vehicle landing footprint must be matched 
with the set of reachable landing coordinates, as deduced from the vehicle 
ranging capabilities exhibited in Figure 3-2. This consideration will then 
define a third range of acceptable encounter conditions. 

The desirable entry path is, therefore, one which is included in the 
three entry corridors defined above. If the intersection of these three 
sets is empty, then preservation of the direct entry mode will require at 
least one of the following 

. modification of vehicle aerodynamic capabilities 

. provision of guidance and control capabilities 

The first item is discussed in the next two sections. The second is considered 
later in Chapter 5. 

3.3.2 Ballistic Coefficient Effects 

The ballistic vehicle aerodynamic characteristics are represented by 
the ballistic coefficient B. Figure 3-8 exhibits the effect (at least on 
down-range) of first order changes in this parameter. However, the strong 
non-linear character of the entry dynamics prevents extrapolation of these 
results to significant changes in B. 

This section considers the impact of significant modifications to this 
parameter, resulting from changes in the vehicle drag coefficient, reference 
area or mass. In particular, we consider a reduction of the ballistic 
coefficient by a factor of four, i.e., a vehicle similar to the baseline but 
with a radius of .5 m, rather than lm. Note that this modification 
implies an attendant reduction in the nose radius utilized for heat 
transfer computations. 

To provide insight into the effect, we exploit the approximate solu- 
tions given in Appendix A. We first consider the entry corridor defini- 
tion. The direct entry overshoot boundary, for a ballistic vehicle, can be 
defined as that value of the encounter flight path angle for which the exit 
velocity is equal to the circular orbital velocity, i.e.. 
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V E = 


( 3 - 10 ) 


The atmospheric braking undershoot boundary flight path angle can then be 
shown to satisfy 


(r e /h $ )(1 - cosy 0 ) = Vg r e /v o ' 1 (3_11) 


and is independent of the ballistic coefficient. Therefore, within the limits 
of applicability of this approximate solution, changes to the vehicle ballistic 
coefficient have no Impact on the entry corridor definition. That this 
solution is not totally sufficient In this region, however, is Indicated 
by the fact that the resultant predicted exit flight path angle Is the 
negative of the entry value which is clearly not the case, as the exit angle 
should tend to zero. 

To assess the Impact of ballistic coefficient changes on peak g-loading 

and heating rates, we can utilize the approximate solution in Appendix A 

for ballistic entry at medium and large flight path angles. The peak g- 

loading (see equation (A-24) does not contain any dependence on the ballistic 

coefficient. The terminal velocity, and therefore the impact shock, however, 

are affected. To preserve the drag-gravity equilibrium, changes in the 

ballistic coefficient are directly offset by changes in terminal velocity, 

2 

such that the product B v^. is constant. 
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The reduced coefficient vehicle terminal velocity is therefore 


v 


f 




(3-12) 


where the superscript b indicates values corresponding to the baseline 
vehicle. 

The peak heating rate is (see equation (A-25)) dependent on the 
ballistic coefficient. A change in this parameter causes a corresponding 
change in the peak heating rate such that 



Thus, in general, the character of the entry path is relatively 
insensitive to changes in the ballistic coefficient. To verify this assump 
tion, numerical results have been generated for the reduced coefficient 
vehicle; the down-ranges are plotted for a shortened range of encounter 
flight path angles in Figure 3-9, and sample values of trajectory parameters 
are listed in Table 3-5 for the max relative velocity entry case. For 
comparison purposes, the down-range curve for the baseline vehicle is also 

included in the figure. 

Comparison of Tables 3-3 and 3-5 shows, as expected, small changes in 
the peak g-load, and approximate doubling of peak heating rate and total 
heat load. Significant increases in the down-range and decreases in the 
time of flight relate to the higher vehicle velocity corresponding to any 
given value of drag deceleration. The down-range deviations, however, 
converge on the baseline vehicle values as the encounter flight path angle 
tends towards the vertical. 

The similar shape of the down-range curves for the baseline and reduced 
ballistic coefficient vehicle imply that similar error sensitivities can be 
expected in the two cases. The major difference between the two curves 
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MAX g LOADING- 







1 - J' 


acceptable landing dispersions can be achieved; i.e., the acceptable entry 
corridor is wider. y 

3.4 Atmospheric Braking 

istirs™ 5 : eC r°" COr,Si,JerS e ' ,a,uat1on of corridor , capture character- 
sties, vehicle stresses and error sensitivities for the atmospheric braking 

maneuver, ,.e„ where the vehicle is passed through the ensphere to 

issipate part of its kinetic energy, yet retains sufficient velocity to 

re-ex,t into a captured elliptical orbit. Furthermore, since the perigee of 

the established orbit is within the atmosphere, the vehicle will (in the 

a sence Of a velocity increment) re-enter in subsequent resolutions, gradually 

dissipating more energy until it undergoes direct entry and lands. The 

gradual nature of the energy dissipation implies lower peak , -loadings and 

method V a tlm CaPtUre m ° de represents a conceptually attractive 

method from a vehicle stress viewpoint. 

Section 3.4.1 considers the path characteristics for a ballistic 
unguided vehicle. Section 3.4.2 assesses the impact of ballistic coefficient 
modifications, and Section 3.4.3 considers the lifting vehicle case. 

3.4,1 Ballistic Vehicle 

He consider here the baseline vehicle, with ballistic coefficient of 
-25 m /<g and zero L/D. 

The entry corridor is limited from below by the condition of vehicle 
re-exit at still hyperbolic velocity, and conseguent loss of the sample, and 
from above, by direct entry to a landing. The entry corridor is then givln 
by ttat range of encounter flight path angles such that the exit velocity 
j u pa ™ 0 ' c ( v ex < 2 u/r at - overshoot boundary), but supercircular 
(V ex > w/r a t ‘ ur| dershoot boundary). Figure 3-12 plots inertial exit 
velocity versus encounter flight path angle for maximum and minimum earth 
relative velocity encounter. The acceptable encounter corridor is thus seen 

exit inenia°i f r T' T ^ W PrMe " ts the “™ pp "«"9 

exit inertial flight path angle. 

The two exit parameters determine the nature of the resulting orbits. 
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The corresponding perigee and apogee altitudes are given in Figure 3-14 as 
a function of the encounter flight path angle for the maximum relative 
velocity encounter case. Also included in the figure is the resulting orbital 
period, as given by 


P = 



(3-14) 


■ * <r p + r,)/2 

Table 3-6 presents sample values of peak g-loading, peak heating rate and 
total heat load encountered within the corridor for the maximum relative 
velocity encounter case. 

Since the vehicle exits the atmosphere in an elliptical orbit with 
perigee within the atmosphere, it will re-enter on subsequent revolutions. 
Furthermore, the inertial re-encounter conditions are fixed by the previous 
exit conditions, i.e., * 


n _ n-1 

v entry ~ v exit 

n n-1 

Y entry ~ ’^exit 


(3-15) 


Subsequent passes through the atmosphere result in either a direct entry to 
landing, or a reshaping of the orbit in terms of a reduction in the perigee 
and apogee altitudes. Figure 3-15 presents the changes in perigee and apogee 
altitudes of subsequent revolutions relative to the parameters of the 
initial orbit, as a function of the initial encounter flight path angle. 

Figure 3-12 exhibits the narrowness of the entry corridor and steepness 
of exic velocity vs. encounter flight path angle curves. It may be expected 
from this, that high error sensitivities will be encountered for this 
maneuver. Figure 3-16 presents the error partial s of exit velocity to errors 
in initial encounter velocity and flight path angle, vehicle ballistic 
coefficient and atmospheric density. Figure 3-17 presents the corresponding 
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Table 3-7: Baseline Ballistic Vehicle Exit Velocity and Flight Path 

Angle Error Sensitivities 
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4.0 Maneuvering Vehicles 


4.1 Introduction 

In this section we consider vehicles with a propulsive maneuvering 
capability. Two cases of application of this capability are analyzed: 

i) establishment of initial capture orbit through thrust decelera- 
tion of the returning vehicle on an hyperbolic orbit. 

ii) Perigee modifications of initial orbit established either 
through thrust deceleration or atmospheric braking. 

The first maneuver can be used to park the returning vehicle in an orbit 
about the earth without incurring the high sensitivities encountered in 
Section 3.4 for atmospheric braking. Furthermore, the resulting orbit is 
stable since the resulting perigee is not forced to be within the atmosphere. 

Following this initial capture, the vehicle's kinetic energy can be 
'•’•ssipated by lowering the perigee into the earth's atmosphere. Proper 
placing of the perigee altitude allows establishment of a controlled rate 
of orbital decay to a set of conditions from which either direct atmospheric 
entry in a more benign environment or orbital recovery by a space transpor- 
tation system (i.e., Space Shuttle) can be effected. 

This second maneuver can also be applied to rectify the orbit resulting 
from an atmospheric braking pass. In this case, however, orbital decay 
control may require either lifting or lowering of the orbit's perigee. 

4.2 Thrust Braking 

We consider thrust impulses applied at the hyperbolic perigee in 
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direction opposite to the velocity vector «t thet point. The impulsive 
velocity chsnpe required to place the vehicle on an elliptic orbit of 
eccentricity e Is given by 


Av - ^2u/r p ♦ - ^ u (l+e)/r p (4-1) 

The apogee radius of the resulting orbit Is then 


r a * r p (1 + e)/0 * e) (4-2) 

Figure 4-1 presents the resulting apogee altitudes above the earth as a 
function of the applied velocity change for two values of perigee altitude: 
h p “ 150 kms and h p - 1000 kms. 

The corresponding orbital period is given by 

P ■ 2* ^a 3 /u 


and Is plotted vs. the apogee altitude for the two values of perigee alti- 
tude In Figure 4-2. 

It Is therefore apparent that modest thrusting maneuvers 
(Av - 600 - 1000 m/sec) are sufficient to capture the vehicle into earth 
orbit. Recovery with space transportation systems, such as the Shuttle, 
however, would require significant further reductions In perigee velocity, 
or, equivalently, significant reductions In orbit apogee altitude. 

^ • 2 ^ Q n Maneuvers 

The additional required deceleration may be achieved, at modest cost, 
by lowering the orbital perigee Into the atmosphere and exploiting the 
resulting drag force. Once the desired vehicle state Is obtained. Shuttle 
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recovery may be effected, or the perigee may be raised to place the vehicle 
mo a stable, but lower energy orbit. Alternatively, the vehicle may be 
allowed to decelerate to the point where encounter with the atmosphere will 
cause it to undergo direct entry and landing. 

A similar scenario can be visualized for the atmospheric braking 
capture mode analyzed in Section 3.4. Following the Initial pass through 
the atmosphere, a perigee rectification thrust maneuver could be applied to 
either raise the perigee out of the atmosphere to yield a stable orbit, or 
simply place the perigee at such an altitude that a controlled amount of 
deceleration will be effected during the next atmospheric pass. 

In either case, application of the velocity change impulse at apogee is 

assumed. For small applied impulses, the resultant change in perigee alti- 
tude is given by 


Ah p = (2P/ir) ^ p /r a Av (4-4) 

The change in perigee altitude due to a unit velocity impulse Av is plotted 

in Figure 4-3 versus the orbit apogee altitude for two values of perigee 
altitude. 


4 . 3 Or b i t a 1 Decay 

Given that the orbit's perigee is within the atmosphere, the vehicle 
will experience a period of deceleration in each revolution. The resulting 
velocity change can be approximately considered as a tangential velocity 
impulse applied at perigee. Its major effect, therefore, will be a reduction 
in apogee altitude, given, analogously to Equation (4-4), by 


Ah a = (2PM ^ Av (4-5) 

In actuality, a reduction in perigee altitude will also result, but 
this effect will be minor except in the immediate vlc^ty of the atmospheric 
braking corridor undershoot boundary. This fact was illustrated in Figure 
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(4-7) 


6Y = -V 2 VV r at 

“here r #t is the atmospheric edge radius. 

The sensitivity of apogee altitude decay to er-o-s in the in *• 

Perigee can then be approximately expressed ' fo- hiohiv l ! 
as J CA ^ rei5ea » to. highly eccentric orbits. 


3(Ah )/9r 

a p 




(4-8) 


Neglecting the variations of ( 3 (av)-/ 3 y ) 
encounter velocity allows using the curves of 


with small changes in 
Figure 3-12 to represent this 
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parameter. Utilizing the additional approximations 


2u 


Vv 


(4-9) 


Equation (4-8) is rewritten as 


3(Ah a ) 

^T 




(4-10) 


The sensitivity 3(Ah a /r a )/3rp is plotted in Figure 4-5 versus the perigee 
altitude. 

The generated sensitivity curve is only truly applicable for approxi- 
mately parabolic orbits. A more complete analysis of these effects would 
require generation of velocity change sensitivity curves of the form of 
Figure 3-12 for the range of possible encounter velocities. This effort 
has not been undertaken in this study. 


56- 


INTERMETRICS INCORPORATED • 701 CONCORD AVENUE • CAMBRIDGE. MASSACHUSETTS 02138 • (617) 661-18*0 


5.0 System Requirements 


Chapters 3 and 4 analyzed three distinct HSR capture and recpvery 
methods : direct entry, atmospheric braking, and thrust deceleration. Each 

of these methods presents requirements on system capabilities in or er 
meet constrained envelopes of desired terminal conditions. 


interactions of the vehicle with the earth's atmosphere were seen to be 
characterized by high error sensitivities, unless substantial vehicle 
stresses, in the form of structural loading and heating rates, could be 
endured. Thus, restrictions on allowable landing dispersions, or on 
acceptable atmospheric re-exit conditions, will generally require the 
provision of a guidance and control capability to maneuver the vehicle 
towards the desired terminal condition envelope. 


on the other hand, vehicle maneuvers involved in the thrust deceleration 
capture mode, or for orbit rectification following the initial atmospheric 
braking pass, imply the provision of an impulsive capability to e ec 
velocity changes required, and of some means to control their direc ion 

and maanitude. 


This chapter analyzes these system requirements. Sections b.l, 5.2. an 
5 3 discuss guidance and control requirements for atmospheric deceleration 
modes, while Section 5.4 discusses thrusting maneuver requirements. 


5 . 1 Guidan ce Bas eline 

Ue are concerned with controlling terminal conditions in the form of 
down-range for the direct entry «de and of exit velocity and flight pa h 
angle for the atmospheric braking maneuver. It should be noted that 
guidance and control capability could be utilized to condition other 
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trajectory parameters (such as vehicle stresses); however, these considera- 
tions are not a subject of this discussion. 

The available control force is assumed to be aerodynamic, i.e., lift. 

The control parameter is then the magnitude of the lift vector component in 
the up direction (e.g., fixed L/D vehicle with roll control of Apollo). 

The two considerations of primary interest in specification of system 
requirements for support of atmospheric interface capture modes, are 

. control envelope 
. mean-square terminal error 

The first consideration assumes perfect knowledge of the state of the vehicle 
and perfect derivation of control commands, and is concerned with specifying 
the control capability, in the form of a minimum vehicle L/D ratio, required 
to i) reach desired terminal conditions from constrained initial states, 
and ii) compensate for uncertainties in initial conditions and disturbances 
along the path. 

Thus, the provision of a control capability may allow the specification 
of a corrmon landing site for the whole family of expected MSR return 
trajectories, which might not bv’ reachable through ballistic direct entries 
due to constraints on return orbit inclination and allowable vehicle 
stresses (see Section 3.3). 

On the other hand, the atmospheric encounter process was seen to be 
characterized by high error sensitivities. Chapter 3 indicated that the 
primary error source in this process is uncertainty as to the atmospheric 
encounter flight path angle. The provision of a control capability thus 
conceptually permits compensation of a range of these errors by steering 
the vehicle to the desired terminal conditions from the offset initial 
state. 

The second consideration establishes the accuracy of the guidance and 
control system required to determine deviations from the nominal path and 
to generate corrective actions. In the presence of uncertainty, it is 
not possible, in general, to meet desired terminal conditions with zero 
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error. Instead, we are forced to consider minimization of terminal errors 
in a statistical or mean-squared sense. 

The certainty-equivalence principle or separation theorem of optimal 
control indicates that the covariance of terminal errors is equal to the 
sum of the estimation error covariance and the control error covariance. 

Given that the covariance matrices are positive semidefinite, it is clear 
that mean-square terminal errors can never be made smaller than the mean- 
square errors in estimating them, regardless of the control logic used. 

We will then assume that the optimal control shall be utilized, and will 
concentrate on the accuracy requirements of the estimator. 

During the atmospheric trajectory, a plasma sheath will form around 
the vehicle, blacking out all external communications. The navigation sensors 
utilized for estimation of trajectory deviations must therefore be self- 
contained to the vehicle. A standard sensor for this application is an inertial 
navigation system (INS). 

We are therefore concerned with terminal condition estimation errors 
of the INS resulting from initial errors at atmospheric encounter (i.e., 
alignment) and instrument (i.e., gyro, accelerometer) errors which disturb 
the propagation of the estimated state along the path. 

5.2 Direct Entry 

We determine the set of reachable states, for given encounter condi- 
tions, by considering the range of terminal conditions between those 
achieved by full -up and full -down control. Figure 3-1 presented resultant 
down-ranges for a ballistic vehicle as a function of encounter flight 
path angle, and Figures 3-3 and 3-4 translated these into achievable landing 
latitudes. Figure 3-10, on the other hand, presented these down-ranges for 
fixed L/D vehicles will full-up lift. A new set of landing latitudes could 
then be similarly obtained. Furthermore, a third family of curves could be 
generated under the assumption of full-down lift, with a corresponding set 
of landing latitudes. The range of landing latitudes, for a given encounter 
flight path angle and orbital inclination, represents the set of reachable 
impact, latitudes for the selected vehicle L/D ratio. 
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From Figures 3-1 and 3-10 it is clear that a broad set of controllable 
errors exists for minimum lift capabilities. For example, assume a nominal 
trajectory with encounter flight path angle of -7.5° at maximum relative 
velocity entry. The resultant down-range is 548 kms. At full-up lift, 
this down-range is achieved with an encounter flight path angle of -9° for 
L/D = .1 and -11.5° for L/D = .2. Therefore, encounter flight path angles 

uncertainties of 1.5° for L/D. .1 and 4° for L/D = .2 are controllable with 
up lift. 


Evaluation of INS terminal condition estimation error sensitivities 
requires selection of a nominal trajectory, since some of the involved 
error sources are dynamic. The selected nominal corresponds to maximum 
relative velocity entry of baseline ballistic vehicle at an encounter fli 

path angle of -7.5°. Characteristics of this entry path were given in 
Table 3-3. 


ght 


The results presented here have been generated utilizing a space-stable 
INS mechanization model, so no sensitivities to gyro torquer scale factors, 
or gyro misalignments will be presented. Furthermore, to obviate the impli- 
cations of relative alignment of instrument axes with the operating drag 
force, the results are presented for general types of error sources, rather 
than for individual components. 


Table 5-1 presents down- and cross-range terminal errors resulting 
from unity INS error sources, including initial alignment errors, 
accelerometer biases, scale factors and misalignments, and gyro biases and 
g-ser-,itive drifts. The presented values thus reflect individual error 
source influence functions onto the INS estimated terminal state. In 
addition, given the short time of flight intervals relative to the Schuler 
period, it can he expected that initial position errors translate directly 

into corresponding terminal errors, and that initial velocity errors integrate 
in similar open loop fashion. 
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Source Influence Functions - Direct Entry Maneuver 






” v 


- short high-deceleration se^ent anVT trajectories consi 

passive, aerodynamical^ controlled vehicle, “(To ^ 

control capability is concentrated in th i < <‘ 9 ** f X * d L/0)a most of the 
" Uh little control force available during the l' ^ d<?Celerat1on se 9">ent, 
alternate view of the control problem is then t aT^ An 

set of equilibrium fall insertion r th * f matchfn 0 « desired 

errors ^cumulating anJZ p „ (nt " S ' “ MCh ,NS «‘<™t,on 

E:H- r — • 

■™ r ‘ ““ " «- 

5 • 3 A tnwsgher j cj |ra k i nq 

The set of reachable state* * fc . 
t;on of Figures 3-12, 3 -l 9 . and 3-13 'l * C ° 

of o,r ves reflecting the ape, 1cat , 0 „ of . ' ! Ge " eratfon of a third family 
determination of the range of exit i ' d °“" -"'c 01 ‘hen allows 

frm a counter flight pat h a „g u ”* *"** M " a " 9 ’“ '—able 

curves, for example/cEsidlrTl 0 "^'' <S e " dent fr01 ” these 

f'*9ht path angle of -3.7" The “T re ’ atiVt ' '' e ’° City —enter at a 

this exit ve.oci^ is « TT' ^ ’* «*" 10 -/«. 
angle of -3.85’ for l/D = , 4 . ftr *" -counter flight path 

conge of controllable errors (for up lift)',- i‘ M ^ L/D ‘ ' 3 ' The 

t/D * •'« ond ,a 5 r f o r L/0 . | 3 0r “ P " H) ’ s — -'5' for I/O . , r for 

nominal trajectory coEl N poL S ^g a to 0 max, CCUraCy haS been CarriM ° ut for a 
the baseline ballistic vehicle at an en'""”" re,4t ' ve velocity encounter of 

Tha ,NS again represents a space "mT '’h^ 1 an9 ' e ° f - 3 ' e5 ”- 

pace st ^o1e mechanization. 


•‘ 63 - 


INTFRMfTRiOS INCORPOR ATPn ?a 

comcoRo*v ENUf . C4MBRinc , f M4SSAChusei|s wim _ (#i7( ^ (<w 


Source Influence Functions at t=200 secs - Direct Entry Maneuver 
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Atmospheric Braking Maneuver 
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6.0 Inertial Subsystem Hechanizations 


In order to perform guided propulsion or aerodynamic maneuvers, it is 
necessary to determine and maintain an estimate of vehicle attitude, and 
possibly specific force, depending on the type of guidance. In this chapter 
some specific hardware alternatives will be described, along with their imple- 
mentation characteristics as regards complexity and reliability. 

6 . 1 Strapdown vs. Gimbal Ter 

6.1.1 Gimballed Inertial Subsystems 

This type of implementation consists of a gimballed stable member, upon 
which are mounted gyro torque sensors and accelerometers. Once the inertial 
orientation has been determined, it is possible to maintain this estimate 
without the explicit intervention of attitude computations or other data 
processing functions; attitude excursions which result in torques applied to the 
stable member are sensed by the gyros, and counter torque commands are fed 
back to motors mounted on the gimbal axes. As the gyro input axes can be used 
to define an inertial coordinate system, output from co-aligned accelerometers 
can be used directly to integrate a state vector maintained in this frame. 


6-1.2 Strapdown Inertial Subsystems 

In this type of system, the input axes of the gyro and accelerometers 
have a fixed orientation with respect to the vehicle body axes. Given an 
initial estimate of body attitude with respect to some Inertial frame, the 
changes In vehicle orientation are incrementally sensed and used to Integrate 
the body-to-lnertial transformation matrix. This matrix Is then used to 
periodically convert accelerometer measurements to the reference navigation 
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coordinate system. Generally speaking the strapdown system offers less 
-hanica, complexity at the cost of additional processing to utilise 

6.2 Gyro Types: Floated vs Dry 
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of these calculations are shown in Table A for direct entry and Table B 
for Atmospheric Braking. 

6.3.1 Di scussion of Error Sources 

Examination of Tables A and B shows that in both cases, the error source 
which dominates the terminal errors is the Initial stable member misalignment. 

This deserves some comment. 

In order to align an inertial system, it is necessary to determine the 
orientation of the stable member with respect to the fixed stars, or equivalently, 
with the coordinate system of the earth plus a knowledge of time. For the 
assumed IMU in the MSR vehicle, this Is accomplished by using a star sensor 
or similar device to establish the direction of two or more stars with respect 
to the vehicle. These directions are then transferred to the stable member 
coordinates system using the glmbal angle readouts, in this case, the one- 
speed resolvers. 
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A resolver Is an electromechanical device and is subject to many random 
non ran om errors. In the case of one-speed resolvers suitable for use 
on a small platform, these error sources will total 2 tQ 3 minutes of 
arc, lo. In addition, the resolver outputs must be converted to digital 
form for use in the navigation computer, thereby creating additional error 
sources, including the quantization error of the resolver to digital converter 
This resolver plus converter error is the principal cause of the initial 
misalignment in our assumed system. 

It should be noted that earth based systems (aircraft, shipborn. etc) are 
not subject to this error source. These systems are normally aligned by 
gyrocompassing, which uses the gyros and accelerometers on the stable member 
o determine the plane which includes both the earth’s gravitation vector 
at the location in question and the earth’s angular rotation vector. The 
me in that plane which is perpendicular to the gravitation vector is taken 
o e t e North-South direction. This alignment is accomplished without 
any need to transfer direction information through the gimbals. 

The problem of in-space alignment of gimballed inertial platforms has 
een attacked many times. In both the Apollo and the Space Shuttle manned 
space systems, multi-speed -esolvers were added to all of the gimbal axes 
This adds considerable complexity to the gimbal structure and to the resolver- 
to- igital converter. Other possibilities also exist, such as: 1) substituting 
a different, more accurate gimbal angle readout. 2) transfer the alignment 
orma ion irectly to the platform optically, bypassing the gimbals, or 3) 
somehow mount the star sensors inside the gimbals on the stable member, 
one of these solutions appear viable, and all would certainly be expensive. 


6,4 jtrapdown Sy stem Navigation Performer* 

With a strapdown system, the gimbals and resolvers are eliminated, along 
with the error contribution from the resolvers. There are other error sources 
however, which contribute to the actual misalignment, including some which are 

not present in gimballed systems. One of these new error sources is the gyro 
torquer scale factor error. 
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In a strapdown system, the gyro must follow the motions of the vehicle. 
This is accomplished by closed loop torquing the gyro to null the error 
output, and sending the gyro torquing information to the digital computer, 
which integrates the incremental torquing information to keep track of the 
attitutde of the vehicle. The difference between the angular rate actually 
produced by the gyro torquing system and the angular rate used by the computer 
to keep track of vehicle attitude is the torquer scale factor error. This 
error contributes to the error in knowledge of the true vehicle attitude as an 
open loop integration, proportional to the total angle through which the 
vehicle rotates. Furthermore, the error terms for torquing of a gyro in 
opposite directions are in general not identical and in some implementations 
are substantially uncorrelated, so that oscillatory rotational motion of the 
vehicle leads to a growing error, not a cancellation of the error. 

A reason*^ e value for the gyro torquer scale factor error in the type 
of system we are discussing is 75 ppm, la. Exactly how this will ca.pl e into 
the errors in initial alignment and navigation performance during entry is 
difficult to predict until a specific operating scenario ij_avai lable, but it 
is worth noting that this error will cause a .47 mr (1.61 mTn) error after a 
2 it radian rotation of the vehicle. 

Other error sources which contribute to initial misalignment of the 
strapdown navigation system are mounting and machining errors between the 
navigation system base and the star sensors, errors in the star sensor itself 
and gyro drift between the time of alignment and arrival at the beginning of 
the entry maneuver. 

In connection with the last mentioned error, it should be noted that 
0.015 deg/hr gyro bias rate assumed for the gimballed system , the alignment 
of the system deteriorates from this cause alone at 0.26 mr/hr (0.9 mTn/hr). 
On the other hand, with a strapdown system, it should be possible in a reason 
able length of time to calibrate the gyro bias during the alignment process 
to a value of .003 deg/hr or better. 
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With gimballed systems, the only practical way to obtain significant 
redundancy is to completely duplicate the platform, inertial components and 
gimbal systems, however many times are required to obtain the desired level 
of redundancy. This is expensive in terms of both weight and power consumption 
especially since each additional system only adds one more level of failure 
tolerance. In other words, if the z gyro fails on one platform and the x 

accelerometer fails on another, both platforms are essentially out of 
operation. 

With a strapdown system, effective redundancy can be added in much 
smaller increments. The stable member (navigation base) of a strapdown system 
is essentially a passive block of machined metal which has an extremely low 
probability of failure. A single, non- redundant, such navigation base can 
be used on the vehicle without significantly effecting the probability of 
mission success and furthermore, this base can easily be made large enough 
to accomodate redundant inertial components. By doing this, the loss of 

one of two redundant z gyro axes and o ,e of two redundant x accelerometer axes 
will not put the system out of operation. 

Another design feature which can effectively increase the redundancy level 
of a system without significantly increasing the weight or power is to mount 
the inertial components with their input axes at skewed angles rather than 
paral.el and perpendicular to each other. By doing this, one component can 
effectively (with some degradation in performance) act as a redundant component 
for more than one other component. This is an area which has received con- 
siderable study, both theoretical and practical, but not all of this work Is 
applicable to the current circumstance since much of it involved use of 
single input axis instruments rather than two input axis instruments. Several 

systems using these principals have been built and flown in space applications 
with good success. 


This ,s an area which should receive additional study directed towards the 
particular current application. This is particularly true because with re- 
dundant hardware the simple value of a Mean Time To Failure is not a proper 
indication of the probability of successful completion of the mission. Instead, 
this probability is a more complicated function of the actual duration of the 
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6.7 Alternatives 


If the navigation performance of the dry gyro configurations described 
in this chapter is found to be inadequate for the intended mission, the next 
step up in inertial component performance would be to switch to the older 
floated gyro technology. The costs associated with this change are 50* to 
1 00* increases in estimated weight, some increase in operating power, and 
possibly a requirement for active temperature control of the gyros when the 
system is net in operation. It should be noted, however, that unless some 
way to increase the accuracy of the initial alignment can be found, there 
is little to be gained by increasing the performance of the inertial 
components. 

6.8 Other Technologies 

Several new gyro technologies are presently being studied and developed 
by a number of investigators, laser gyros and atomic gyros are included in 
this group. While these new technologies appear to have some promise, none 
of them are currently sufficiently mature to be baselined for the present 

application, and it appears that they will not be mature for a least another 
3 to 5 years. 

6.9 Limitations on Mission Profile 

There are several practical limitations imposed on the mission profile 
by the requirements of the inertial instrumentation described in this chapter. 
Among these limitations are: 

1) in order to obtain the assumed accelerometer performance, the peak 
acceleration during tho mission should not exceed 10 to 15 G. 

2 ) If the strapdown system is to be used, the limitations in gyro torquer 
capability will restrict the maximum angular rates of the vehicle to 
200 to 400 degrees per second, and average rates over one second or 
more to 100 degress per second. 

3) If the strapdown system is to be used, spin stabilization of the 
vehicle should be avoided If possible, If spin stabilization must 
be used, the strapdown design may not be applicable. 
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7,0 Risk Analysis 


A quantitative analysis of the reliability of various recovery methods 
is possible in principle if a complete description of the critical systems is 
available, together with the failure modes of each component, the probability 
thereof, and a statistical description of the characteristic output of each 
component, given that it has failed. In fact, such information is generally 
not available with respect to similar equipment used for commercial purposes 
(such as inertial measurement units), and in any case, the rate of failure is so 
small that insufficient quantitative data are available to characterize the 
failed outputs. This chapter will therefore not attempt a detailed quantita- 
tive estimate of the reliability of various methods, but instead identify the 
key system components necessary for successful performance, and the effect of 
various failures or each of these. Such effects depend not only on the nature 

° * e fanure ' but the fu 9 ht regime of the vehicle at the point where the 
ailure occurs. Possible outcomes of a given earth encounter sequence are taken to be 

a) Recovery of the sample module according to the constraints and 
protocol laid down for preservation of physical and biological 
integrity of the container. 

b) Sample recovery involving violation of such constraints and protocol. 

c) Loss of the sample module in the biosphere. 

d) Loss of the sample module outside the biosphere. 
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the integrity of the entry module relative^toth 9 ^ !*"“ teSt1 " 9 c,n " prove " 

^ will be discussed in th ' ' 

tenure mode connected with ,m LT?? 9 ^' *•" re ™'" s a " “"“voidable 
fatlure Of the bus command andcontro, la^Mlit fr “" 7 bus - or 

entry. This mode is common to all rern re!uU, " s f " “"controlled 

"ooule to enter or -,ss through the ter " MCh lhe 

rt fS bb » p ved ‘ha, the unguided entry veht atm ° SPhere - ‘ b <e respect 
tt limits the period of exposure to « , '* preferab,e to all other, as 

deployment, rather than all the way ’ Prior to 

the way through the aerodynamic phase. 

7 ' J ^aested _0ptimum Basel,-.,. 

The analytical results presented i, fh , . 

recovery modes based on technical infeaslbilit" p"™' ‘" y ° f tbe 

JSC and JPL personnel concerned with t.J ' ’ ^ dt “ussions with 
>t seems clear that the preferred method f P ' yS ’ Ca ! > '" 1 bio, °5 (c al integrity, 
technically feasible and ofcers ma i ° reCOVery “ ,n be that which is 
Pf tha terrestrial environment. F«.“ '"T" aSatnSt baplt - pp "tamination 
faded that the recovery method should not i„ r"”" 1 Prev,0 “ s ' y * ft ,s 
module to enter or pass through the at k V ° V6 targetln9 the bus or sample 
a ‘“dy Will show that an earth orbita- c t " “ ^ tba ‘ P “fher 

tan be designed such that there is no mea P s' ,re hl teCh '' f< " ,e ’ US '" 9 retropro b“'*fon, 

“tth the biosphere fur any failure 7 ' prpbab "'‘y of encounter 

eny failure mode, including uncontrolled burn of the 
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Sj,stem re,u,>ewnts to ^ ««* 

” ***»«* “termination can be by star tract, t ' ind Wrf °™"« 
,n use for Similar missions, P ch "°'° 9J ' cur rently 

rr: 

conditions 'and^gn'tiorTtime' * C °" ,r °" ed * C0Bt " , ""» app ™ p * 

Operational advantages include the permissibilitv of , 
orbital insertion maneuver even when the , , * POrformlng the 
Integrity is tnown to be ”’ P * Ct #P M °’°’ ! ‘ a ’ 

and trajectory control requirements- and th • - ° rbU det ormi:iation 

or On-orbit Investigation of the sa ' * poss,tlf ’tty of delayed recovery 

of the sample to a terrestrial laborV *° dUl ' S " d C “ n,ents - Transfer 
recovery to Shuttle operational ,»■! 'T > "" ** accompl, 'sbed by remote 
* Placed a Shuttle -bor’ is 1 ’ ^ Sa " P " TOdp ’ a p an 

-tastrophic outcome of a Shuttle entry'”' *° “ ithSta " d ^ 

sequence win be discussed in the next chapter aSS ° C ' ated " nh tb ' S 
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8.0 Orbital Recovery Operations 


Thir chapter discusses the operational requirements to effect recovery 
of a sample module to a Shuttle payload environemnt. The most reliable flight 
vehicles available today are commercial passenger aircraft- These incorporate 
redundant mechanical and avionics components, and are designed to withstand 
stresses well above those to be encountered in normal conditions. As an aero- 
dynamic vehicle, the Space Shuttle is subject to the same flight hazards as 
any manned aircraft, and to the additional hazard that it is a statically 
unstable vehicle through most of its flight regime, and unlike current airliners 
cannot be controlled without continuous computer assistance. As it is an 
unpowered vehicle, its navigational aids are flight critical except in the 
immediate vicinity of its landing field, under VFR conditions. The point 
of these observations is not to disparage Shuttle reliability, tut to a^gue that 
this vehicle is not obviously more reliable than a commercial cargo aircraft. 

It therefore seems likely that special provisions for housing and protecting 
the sample module as a "hazardous cargo" .'ill be required, just as they are 
with such materials in the case of domestic carrier. 


8.1 Pi rect Recovery 

In the case where the earth orbit injection module has further maneuvering 
capability, it may be possible to lower the apogee to within range of the 
Space Shuttle. This is currently expected to be in the vicinity of 
(R + R p )/2 i 225 n.m. without additional OAMS fuel in the payload bay, or 
about 300 n.m. with. In order to reach these altitutdes, the ground track of 
the bus must pass almost due east through KSC. The Shuttle currently provides 
rendezvous radar tracking on non-cooperating vehicles out to about 30 n.m., 
and to about 400 n.m. with a transponder. Although the transponder design is 
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available, there is currently no plan to develop one, and its provision 
is therefore a responsibility of the user. Use of the skin-tracking mode only 
is undoubtedly feasible with the Shuttle on-orbit navigation capability expected 
to be available at that time, provided that the target can be accurately located 
by ground tracking. The remaining question involves actual acquisition of the 
sample bus and placement in the payload cannister. This cannot be done with 
existing technology unless the bus is attitude stabilized and not spinning. 

The problems of achieving such stability or devising mechanisms for recovery 
of a non-cooperating bus, need further study. 

8.2 Remote Recovery 

At the present time, the best forseeable candidate for a remote recovery 
system is the space tug, not yet under development. If this is not 
available, a modified Interim Upper Stage is the alternative. 

The Space Shuttle provides capability to perform rendezvous guidance and 
navigation, braking and station-keeping, as well as the payload manipulator 
system. All of these functions need to be incorporated in any remote recovery 
system, and the current plans for the tug or the IUS do not support this level 
of capability. It is the opinion of the authors that systems for performing 
remote retrieval will become available within the next ten years as an expan- 
sion of the overall capability of the space transportation system. If the 
remote retrieval option is decided upon, an immediate analysis of the per- 
formance requirements for the MSR mission should be undertaken, the results 
of which should be used to lend impetus and direction to the development 
of remote earth orbital operational capability. Bearing in mind that the 
cost, of development can be spread over a number of customers, and the fact 
that remote retrieval is required in the case where a system failure eliminate 
the direct recovery option, this would seem to be the method of choice. 
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9.0 Open Issues and Future Studies 


The most open issue associated with MSR earth recovery strategies is, 
and will continue to be, the level of reliability against back contamination 
which is necessary to gain acceptance of the mission objective. This report 
has argued that on the basis of system simplicity, performance requirements 
and weight, a direct entry, high flight path angle vehicle is the method 
of choice for a minimum system. This report has not addressed the question 
of structural, and hence weight, requirements for a vehicle which must pass 
through a 300 peak g profile, and impact at a speed of greater than 20 fps. 
Furthermore, the question of whether the earth encounter and deployment phase 
can be designed so as to offer an acceptable probability of uncontrolled 
entry remains unresolved. It is therefore recommended that the following 
two areas must be further evaluated for the minimum MSR: 

a) Estimate of the weight of structure and thermal protection system 
for a high angle entry vehicle; and 

b) Preliminary design of the earth encounter trajectory and deployment 
mechanism to minimize exposure to pre-deployment systems failures. 

With respect to the recommended optimum baseline, the two most pressing 
issues are the conjecture that a guaranteed safe encounter/retrofire sequence 
can be constructed, and that suitable mechanisms will exist to recover a 
sample module above Shuttle operational altitudes. These areas can be further 

elucidated by 

(i) A preliminary trajectory analysis of the approach path and retro 
sequence, to show no on-board failure can result in encounter with 
the earth atmosphere; and 
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( 1 1 ) A survey Of potential needs for, and users of, remote operations 

and equipment for earth orbital applications. This survey should be 
oriented toward identifying common functional and performance 
requirements, and should result in a government/industry development 
effort. 
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APPENDIX A 

mathematical problem statement 

A -1 Equations of Motion 

Assuming that the vehicle possesses a ctaKio ™ 
tl,en all operating forces M . ! ,.! 6 C0 " f ’ 9urat1 °" at zero roll, 

contained in the plane defined ^r^C^rT 9 * 0 *" 

to eartn-center vectors. TOe action of the vehicle re, t^ ZTZZT'' 
et;rth is described by [ 1 ] ipnericai 


A-l 
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“ = - \ pv 2 CpA/m - g siny 


^ = \{~r " g)cosy + | pv C,_A/m 


— = v sin y 
dt 


IB = _ e v cos y 
dt r 


P = p(h) 


(A-l) 


where 


v 

P 

y 

h 

r € 

r 

C [ 

A 

m 


CpA/in 


;p v 2 B 


= earth relative velocity 
= atmospheric density 

= flight path angle between velocity and local horizontal 

= altitude above the earth 

= earth radius 

= vehicle radius - r g + h 

= vehicle drag coefficient 

= vehicle reference area 

= vehicle mass 

= B = ballistic coefficient 

-• vehicle lift coefficient 

= D = drag acceleration 


1 2 

-pv C. A/m = L = lift accelerati 
2 L 


on 
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L/D = lift to drag ratio 

R = range = distance travelled over earth surface 

Given the vehicle capabilities, as expressed in the ballistic coefficient 
B and the lift-to-drag ratio I/O, for example, the characteristics of the 
atmospheric environment, as expressed in a functional relationship of 
atmospheric density p to altitude h, and appropriate initial conditions for 
v, y» h, these equations can be integrated in time to define the vehicle 
state at any desired instant. 

The structural loading on the vehicle can be readily computed as the 
sum of the operating lift and drag forces. Its magnitude is simply 


F = moVl + (L/D) 2 (A-2) 

The stagnation point convective and radiative heat rates are given 
by 

~ = k 1 (R N ) -1/2 (p/p o ) 1/2 (v/1000) 3 - 15 

(A- 3) 

~ = k 2 (R N )(p/p o ) 3/2 (v/10000) 12 ’ 5 


where 

Q„ * convective heat flow 
c 

Q r = radiative heat flow 

R N = nose radius 

v = velocity 

p = atmospheric density 

p Q * reference (surface) density 
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Except for extremely high velocities, the convective heat transfer Is the 
dominant effect. 

The total heat flow over the path Is then obtained by Integrating the 
heat rates In time and summing their contributions. 

A. 2 Initial Conditions 

Integration of Equations (A-1) requires specification of Initial condi- 
tions for the vehicle's earth relative velocity v, flight path angle y, and 
altitude h. 

The altitude at which aerodynamic flight Is Initiated Is defined to be 
that of the outer edge of the earth's sensible atmosphere. This boundary 
Is selected to be at an altitude above the earth's surface of 121.92 kms 
(400,000 ft), at which point the atmospheric density Is of the order of 
1 .87x10“® kgs/m 3 . 

For the Mars return hyperbolic orbits, the vehicle velocity, (relative 
to a non-rotating earth), is given, in terms of the orbit hyperbolic excess 
velocity v^, by 

2 1/2 

v^ * [2w/r ♦ O (A-4) 


where 

14 3 2 

p * 3.986012 x 10 m /sec 8 earth gravitational constant 

The corresponding flight path angle y^ Is most easily expressed In terms of 
the aim radius, r of the orbit relative to the target planet 

fl 

tan Yj * ^ (v, 2 r 2 /v. 2 r # 2 ) - 1 (A-5) 

The geometry Is Illustrated In Figure A-1. The aim radius, which Is 
equal to the semi-minor axis of the approach hyperbola. Is given in terms of 
the desired vacuum perigee altitude, rp 
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Equations ( *-1) involve, however, atmospheric airmass relative velocity 
and associated flight path angle, rather than inertial quantities. As 


* * *1 - 2*1 (A-7) 

where ills the earth's angular velocity, the parameter sets are related by 

| 

V * v i V* “ 2(nr/v i )cosy i cosLs1n(|» f ♦ (flr/v^ 2 cos 2 L 
slny * (v^/v)s1ny. 

with 

L * latitude 

| ^ ■ Inertial velocity azimuth from north 

i 

I The possible combinations of entry latitude and azimuth, however, are 

constrained by the relation 


sln^ « -cos 1/cos L ( A _ g 

where 1 Is the angle of Inclination of the orbit to the earth's equatorial 
plane. Noting, furthermore, that, for hyperbolic orbits, v, » q- then 
for all but the steepest values of Y r 1 1 


v " v fO + (ftr/Vj)co$Yj cos 1 ) 


(A-10) 


The Inclination of the orbit Is limited from below by the declination 
of the hyperbolic asymptote, so that the specification of a value of cos 1. 


A- 6 


hi 


E RICS INCORPORATED • 701 CONCORD AVENUE ♦ CAMBRIDGE. MASSACHUSETTS 02138 • (617) 661-lfi 


tit ' ulo md 


and therefore v, is not completely free. 


Sample Mars return orbits, as described by v^, are given in Table A-l. 
Table A-l: MSR Sample Orbits 



v, (km/secT 

x(deg) 

S(deg) 

1 

3.07 

200.02 

27.13 

2 

3.24 

232.58 

-52.45 

3 

2.82 

314.52 

-43.44 


where X is the right ascension and 6 the declination. The corresponding 
atmospheric encounter velocities are computed from Equation ( A-4) and 
Equation (A- 10), and given in Table A- 2. 

Table A-2: MSR Atmospheric Encounter Velocities 



v^(kms/sec) 

W (k " s/sec) * 

^nllons/sec)* 

1 

11.49216 

11.7012 

11.27927 

2 

11.53876 

11.6823 

11.39341 

3 

11.42793 

11.5987 

11.25453 


*Note: for y-j * 0 
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A. 3 Approximate Analytical Solutions 


A number of approximate analytical solutions for the atmospheric entry 
dynamics have been formulated by various Investigators, and have been treated 
in unified fashion by Loh [3], First, the relationship between atmospheric 
density and altitude is written as an exponential 

-h/h 

P s P s e (A-ll) 

where 


p $ * reference (surface) density 
h $ = atmospheric scale height 

Then, noting that 


dt = ' (p/h s ) ft (*-’2) 

and manipulating Equations (A-l), the following pair of equations is 
obtained. 


d COS Y 

'.dp 


1/LuVu . 

2 ( D )( ir )h s + 




d(v 2 /qr) 

dp 


fN . ,V 

rp m 


) 


lv^ai) h 

sin s 


(A-13) 


It is then found that 
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i) The term cosy(l - gr/v £ )h s /rp is relatively insensitive to 
variations in y and p along the path. 

ii) 2/rp « (CpA/m) (v 2 /gr^iny 
Equations (3-13) can then be approximately solved to yield 

COSY - COSy o = N(p-p o ) 

£n[v 2 /v Q 2 ] x (C D A/m)h s (y - y Q )/N (A-14) 

N = L/D) (CpA/m)h s + cosy(l - gr/v 2 )h s /rp 


where p Q , v Q , y Q are initial conditions. These equations can then be solved 
to yield two of the variables (e.g., v, y) as a function of the third. 

Furthermore, the above equations can be utilized to obtain the max 
deceleration and max heating rates for the trajectory. The max deceleration 
point occurs when 


d /dvv .. d , 2, 

d-v ( dt> ' S (pv 1 = 


(A- 15) 


Utilizing Equations (A-13), the maximum deceleration occurs approximately 
when 


(CgA/m)h /p = siny 


(A-16) 


This additional relation allows complete solution of Equations (A-14 ) 
for p, v, 6 at the max deceleration point, and by back substitution, the 
corresponding max deceleration value. 

Similarly, the max convective heat rate is obtained when 


dA, .. d , 

dv ( dt } - d“v (p 


1/2 v 3.15 } 


(4-17) 
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A similar procedure as above allows determination of the max heating 
rate, and the values of the trajectory parameters when it occurs. 

A * 31 illi cit Approximate Solutions - Direct Entry . Further simplifica- 
tions, reflecting constraints on the character of the entry path, allow an 

explicit analytical solution to be obtained. Two examples are provided 
below: 

A. Gliding Entry at Small Flight Path Angles 

Assuming cosy *cosy q * 1, Equations (A-14) can be simplified to 

v 2 /gr e * 1/[1 + (r e /2)(L/D)(C D A/m)p] 

siny * (2h s /r e )/(L/D)(v 2 /gr e ) U ‘ 18) 

The range, or distance travelled along the path, which, for small 
angles of inclination, is equal to the horizontal range, is 
obtained by integrating 

R = f dR = f vdt 

2 (A ~ 19) 

R - (r e /2)( L/ D)£ n C(T-vVgr e )/(l-v o 2 /gr e )3 

The vehicle deceleration is given approximately by 

dt * " D * " g(1 ' v2 /9 r e )/(L/D) (A-20) 


For this trajectory, no peak deceleration is encountered; rather, 
the deceleration asymptotically approaches the value 


-(j?) - g/(L/D) 

Ql max 


(A-21 ) 
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(A -22) 


The peak convective heat rate Is given by 

dQ 7J7' r 

^ 8/R « 8(L/D » <’ r e> 

b * — 11ist * c Entry at M edian and Large Flight Path Anglos 
In this case, the entry path is approximately linear, with 


cosy « cosy 


-(CpA/nOi^fp-PoJ/sinY 


(A-23) 


2 2 
v * v 0 e 


The max deceleration is approximately given by 

( 

* -(1/2 e h )v 2 siny 
max 50 0 

and the max convective heating rate by 


(A-24 ) 



max 



siny o /3e B h $ 



(A-25) 


A - 3 ’ 2 ^XRlic it Approximate Solution s - Atmospheric Braking . For relatively 
small variations in y and p throughout the path, the parameter N in 
Equations (A-14) can be assumed to be approximately constant. A good repre- 
sentative value is stated by Loh to be that which obtains at the minimum 
path altitude (i.e., y * 0). 


C0SY * ' ♦ N m<» ' Pj 

,1,,! _-(CpA/j|)h s Y/N K 
m 


(A-26 ) 
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where the m subscript indicates values at the minimum path altitude. In 
these equations, negative flight path angles pertain to the path prior to 
the y = 0 point, with positive values thereafter. 

Furthermore, from this representation of the path, the atmospheric exit 
conditions can be derived as 


cos y E * 1 * N m °m => \ z ’ y o 
2 2 (CpA/m)h s Y 0 /Nm 


(A-271 


v r * v m e 
t m 


-> V E “ % /v o 
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